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1 Introduction 



Current experiments have given us no evidence for physics beyond the Standard ModeL 
Indeed, the Standard Model may be vahd, as an effective theory, up to very high energies, 
say the Planck mass, 10^^ GeV, if the Higgs mass is below several hundreds of GeV. Even 
this case, however, does not necessarily imply that no new phenomena will be seen in the 
multi-TeV range which will be explored by future colliders such as the CERN LHC. There 
is the intriguing possibility that, above a parton-parton centre of mass threshold in the 
multi-TeV range, the cross section for the nonperturbative production of many, 0{a^), 
weak bosons W(Z) may be observably large. Unfortunately, there is only circumstantial 
evidence for this to happen which is, to a large extent, just based on the observation 
that leading order perturbative calculations for the production of 0{a^) weak bosons 
violate unitarity near the threshold of 0{a^mw) TeV. This happens both for processes 
with and without baryon and lepton number violation (BLNV). At present 

it is an open question whether the actual (beyond perturbation theory) multi-W(Z) cross 
sections become observably large at such multiplicities and energies (for recent reviews 
see 0-3). New theoretical methods are needed to answer this important question. 



In the meantime one can contemplate the prospects of settling the issue experimentally. 
Multi-W(Z) processes at the LHC would be clearly distinguishable from any other Standard 



Model process [|T0|]- This is due to the hadronic and leptonic decays of the W's and Z's, 
which lead to hundreds of charged hadrons and photons with transverse momenta in the 
GeV range, and to tens of prompt leptons with transverse momenta in the tens of GeV 
range. However, the question whether BLNV can be demonstrated in such a multi-particle 
environment has not been answered conclusively |TD|,|rT[]. It is the purpose of this paper to 
investigate this question. 

Using the results of leading-order BLNV calculations, it is possible to create a phe- 
nomenological model of these processes, with variable parameters to represent the main 
sources of theoretical uncertainty. In this paper such a model is postulated, and imple- 
mented using a Monte Carlo event generator. The program, which takes the form of a 



package HERBVI ||I2| operating within the environment of the HERWIG event generator, 
also allows the generation of B and L conserving multi-W(Z) events. With this pack- 
age computer simulations can be performed which permit the feasibility of demonstrating 
BLNV to be investigated. 

The layout of this paper is as follows. In the remainder of this section we describe the 
theoretical basis of electroweak BLNV. Section]^ is a discussion of our approach to modelling 
BLNV processes at supercolliders, based on the (rather uncertain) theoretical expectations, 
and of practical issues that arise in the construction of an event generator. We formulate 
two contrasting BLNV models, one based on leading-order matrix elements (LOME) and 
one using a simpler phase-space approach with fixed boson multiplicity. We also discuss 
the possible backgrounds to BLNV, and how we model the dominant background, which is 
expected to consist of B and L conserving multi-W(Z) events. 

Section ^ contains the results of our Monte Carlo studies. We present results for the 
simpler model at three energies, including that proposed for the LHC. The LOME model 
requires a higher energy and so in that case we show results only at 40 TeV. We present our 
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conclusions in Sect. ^. Brief details of the HERB VI package, used to obtain the majority 
of the results in Sect. ^, are contained in an appendix. 



1.1 Instanton-induced B and L violation 



Owing to the chiral anomaly |T^-|T5| and the V-A structure of weak interactions, B and 



L are not strictly conserved in the Standard Model |]16|,|T^]. In the presence of nontrivial 



SU(2) gauge fields Wi, the fermionic quantum numbers change according to 

ALe = AL^ = AL^ = ^AB = -AiVcs, (1) 

where 

Ncs^^J d'x e^^' tr (^F,,W, - ^W^.W,^ (2) 

denotes the Chern-Simons number of the gauge field, Fij is the SU(2) field strength, and 
aw = g"^/ (47r) the weak fine structure constant. As is suggested by eqs. (P and @, strong, 
nonperturbative gauge fields, of 0{g~^), are needed in order to change the Chern-Simons 
number, or, equivalently, the fermion numbers, by an integer amount. This is reflected by 
the fact that there exists, on topological grounds, an energy barrier |T^,|T5 between gauge 



fields whose Chern-Simons numbers differ by an integer. The minimum barrier height is 
given by the energy of a static saddle-point solution, the so called "sphaleron" ppj, PT[], 
which depends slightly on the Higgs mass, 

Msp = 2 B{mH/mw) ^ ^ 7-14 TeV , (3) 

aw 

where the parameter B is restricted to lie in the range 1.57 < i? < 2.72. At low energies (-C 
Mgp), therefore, anomalous BLNV processes are only possible by quantum tunnelling under 
the topological barrier, i.e. the corresponding amplitudes are exponentially suppressed by 
a tunnelling factor, 

^i<<M.. oc e-^-/-- ~ 10-^«, (4) 



which leads to unobservably small cross sections or decay rates [|16|,[T^ 



The amplitude is expected to be enhanced when the process involves a large number 
of gauge and/or Higgs bosons. Consider the following BLNV process, which might be 
observable in high-energy proton-proton collisions: 

q + q^7q + 3l + nB W(Z) + H. (5) 



The amplitude for this process can be estimated by means of the instanton approach |ll6| , p!7[] . 
The instanton is a time-dependent tunnelling solution which interpolates between gauge 
fields whose Chern-Simons numbers differ by one unit and passes through sphaleron-like 
fields for particular values of its collective coordinates. It was found in refs. [0,0, by 
expanding the path integral about the instanton, that, in leading-order (LO) in the coupling 
aw, the amplitudes for the processes grow with multiplicity and parton-parton centre 
of mass energy like 



W"^^""^ (6) 
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At high energies, these pointhke S-wave amphtudes violate unitarity. Owing to the facto- 
rial growth of the amplitudes (|^) with the number of produced bosons, this violation of 
unitarity sets in at multiplicities of 0{a^) and at parton-parton centre of mass energies 
of 0{a^mw), i.e. the sphaleron energy (^). 

The corresponding parton level cross section, ctlo, due to the leading order matrix 
element (LOME) of instanton-induced BLNV (^, is given by 
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(7) 



In this expression is the number of vector bosons and n 
coupling constant Q is 

g = 1.6x IQ-^^^ GeV-^^ 



n-B + nn- The effective 

(8) 

V = 2mw/g is the Higgs field vacuum expectation value, and C* is a numerical factor 
representing the effect of averaging and summing over incoming and outgoing states of 
the fermions. The leading order result ^.o; taken at the dominant multiplicities, is 

plotted in Fig. |1[ As anticipated, it rises rapidly with centre of mass energy, and violates 
the S-wave unitarity bound (shown as a dashed line on the same plot) 
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at an energy of order the sphaleron energy. 



(9) 



A violation of unitarity is, of course, unacceptable and indicates the importance of 
higher order corrections. There are strong arguments that the higher order corrections 
to the fixed multiplicity amplitudes exponentiate in the total cross section of BLNV, such 



that, to exponential accuracy, the latter can be written as P^-^T 

nB,njj L ^ 

where 



aw \Mq 



Mo ^ V6n^ 
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(11) 



is of the order of the sphaleron mass (|^). The exponent itself is known only in a low energy 
expansion whose first few terms are given by p^,p8|-|37[l 



F(e) 
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where e = v^/Mq. From this result the following conclusions can be drawn: (i) The total 
cross section for BLNV, though small, is growing exponentially at {mw <^)V^ <^ Mq; in 
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Figure 1: Quark- level cross section for leading order BLNV process, solid line. The S-wave 
unitarity limit is shown dotted line. 



this energy region, the cross section is insensitive to the Higgs mass, (ii) The different 
terms in the perturbative expansion of the 'holy grail function' F become comparable in 
size, and the perturbative expansion breaks down, at \/5 ~ Mq. 

This breakdown of perturbation theory is not surprising. BLNV processes involve the 
electroweak gauge fields penetrating the topological barrier between different vacuum sec- 
tors. As the energy approaches Mq, the height of this barrier, the possibility of crossing the 
barrier classically arises, and the expansion based upon the instanton tunnelling solution 
ceases to be valid. Unfortunately, nothing is known about the behaviour of the holy grail 
function around or above the sphaleron scale. Unitarity arguments [^^^1 suggest that its 



increase will stop at values of order F ~ —0.5, leading to unobservably small cross sections 
for BLNV. However, this question is not settled finally. 

To end this section, we note that at low energies the total cross section (|10D is dominated 
by multi-W(Z) production (n^ ~ otw) rather than by multi-Higgs production. In particular 
one finds at low energies 

SB~^rv+o(.^)). (13) 

47r 3 



3avi/ \8 

whereas 



nn ~ — T^^e' • (14) 
aw 3/ 



1.2 Sphaleron decay 

It is interesting to compare the above results for the instanton-induced process (|D with 
those from sphaleron decay. Since the sphaleron lies atop the energy barrier through which 
tunnelling occurs, it is expected that the lowest energy classical field trajectories which 
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violate B and L will pass (in some sense) close to the sphaleron configuration. Thus, we 
expect the behaviour of the sphaleron to give us some indication of what should be observed 
in an actual BLNV event. In particular, the decay products of the sphaleron (which is by 
nature unstable) should be similar to the gauge and Higgs particle byproducts of a BLNV 
event. Therefore, one again expects [^,0 that the dominant BLNV processes will involve 
0{a^) weak bosons, simply because sphaleron-like intermediate states will typically decay 
into many W's and Z's. 



The decay of the Klinkhamer-Manton sphaleron has been studied numerically in refs. E3 



44|. The sphaleron field configuration was first discretised, and then allowed to evolve under 



the classical Euler-Lagrange field equations with an imposed condition of spherical sym- 



metry |]45|- ^] . The restriction to spherical symmetry is permissible because the sole decay 



channel of the Klinkhamer-Manton sphaleron is spherically symmetric it should be 
noted, however, that this approach may be used only in the limit of vanishing weak mixing 
angle 6'w As a consequence of the imposed spherical symmetry, it is expected that equal 
numbers of W~^, W~ and Z° of any given polarisation state will be produced after the decay. 
This is due to the fact that the spherically symmetric ansatz mixes spacetime and group 
indices, so that particle identities may be interchanged by appropriate spatial rotations. 

As a sphaleron decays and the energy carried by the gauge fields moves away from 
the sphaleron's initial position, the amplitudes of the gauge fields show the 1/r attenuation 
characteristic of spherical radiation. Some time after the sphaleron decay, the amplitudes of 
the gauge fields become small enough that the nonlinear terms in the gauge field equations 
of motion are negligible. Thus, a long time after the decay, the various components of the 
gauge fields behave approximately like free massive fields, and they evolve independently 
of one another. The multiplicities can then be measured by interpreting the classical fields 
in terms of coherent states PD|-P|. 

Figure || shows the expected particle multiplicities in sphaleron decays for several values 
of mlj/m^r. The triangles indicate the expected numbers of transversely polarised weak 
bosons per helicity state and per particle type (n^^), and the squares show the expected 
numbers of longitudinally polarised bosons per particle type (^l)- The crosses show the 
expected numbers of Higgs bosons {nu)- The total number of gauge bosons is then 

TLB = Gtlt + 3fii . (15) 

The expected number of gauge bosons is seen to vary weakly with the Higgs mass and is 
of 0{a^), as anticipated. 

Setting = M^p leads to 

e.p = ^ = -^B{mH/mw) ^ 0.4 - 0.7 , (16) 



which may be computed readily from the table of M^p versus rriH provided in |21]. Figure ^ 



compares the expected gauge boson multiplicities computed from (|T3p, shown by triangles, 
with those given by Eq. (0), denoted by diamonds. It can be seen that the multiplicities 
predicted by ([T3|) are roughly 30% smaller than those expected from sphaleron decays. This 
discrepancy is not very surprising in view of the fact that the perturbative expansion of 
F(e) does not converge rapidly for values of e so near to unity. 
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Figure 2: Sphaleron decay products as a function of the Higgs mass. Points shown on this 
graph show the numbers of each kind of particle in each possible polarisation state: triangles 
represent transversely polarised gauge bosons, squares show longitudinally polarised bosons, 
and crosses denote Higgs particles. 
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Figure 3: Mean number of gauge bosons, n^, expected to be produced in a baryon number- 
violating event, as a function of e = V^/Mq, where Mq = \/QTTmw/aw- Triangles show n-Q 
from sphaleron decays, and diamonds show from Eq. ([T3|) . 

2 Modelling B and L violation 



In the absence of a more complete theory at energies of order the sphaleron energy, we will 
use the gross features of instanton and sphaleron calculations, such as the threshold-like 
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behaviour of the cross section and the parametrically large vector boson multiphcity, as 
a guide for phenomenological modelhng of electroweak BLNV. The steeply rising cross 
section means that BLNV processes have a threshold nature. Three parameters are used 
which encompass the main features of the processes: a threshold energy \/so, a threshold 
cross section ar, and the boson multiplicity, ub- The parton level cross section is then 
modelled as a step function 



a (Vsj =aTe i^Vs - ^soj (17) 
following the approach used in [^. We assume that in the rest frame of the interacting 



partons the events are isotropic, and that the outgoing particles have an energy distribution 
that fills the available phase space uniformly. We first discuss the model parameters and 
general assumptions, and then, in separate subsections, explain the detailed modelling of 
events and background. 

The magnitude of (Tt is obviously a crucial quantity, as it determines the rate at which 
the interactions will occur. It is also difficult to estimate because of the theoretical un- 
certainties, discussed above. For this reason we have performed estimates in terms of the 
number of events required to verify the existence of BLNV, assuming a comparable back- 
ground of B and L conserving 'multi-W' events with a similar boson multiplicity. The 
rationale behind this model for the background is presented in Sect. 



The high energy processes under consideration mean that, for any pp collider proposed 
for the near future, the threshold energy y/so is a significant fraction of the total beam 
energies. Consequently, events will be at large momentum fractions x, and also will tend to 
be central within the detector. The exact size of a/Iq is not significant for our studies. We 
have set this parameter at the point where the instanton calculation violates the S-wave 
unitarity bound, ~ Mq, except for the studies performed at 17 TeV where the value of 
5 TeV (~ 2'mw/c(w) is used. 

Convolution of ( [T7| ) with the appropriate proton structure functions allows a comparison 
of the relative rates of processes with similar parton-level cross sections mediated by the 
possible initial state combinations of qq, qq, and qq at a pp collider. These relative rates 
are plotted versus threshold energy in Fig. ^. The results shown were obtained using the 
EHLQ structure functions for the proton |5^, but the dominance of qq processes at high 



threshold energies is not sensitive to the choice of structure functions. It follows simply 
from the fact that the valence u and d quark distributions dominate over the sea quarks 
and gluons at high values of the parton momentum fraction x. 

For BLNV, which can only occur for incoming qq or qq, the latter contribution is two 
to three orders of magnitude smaller than the former. Therefore the contribution of BLNV 
with the opposite sign (creation of baryons as opposed to antibaryons) is negligible. 

In this paper, two alternative methods for generating the boson multiplicity ub are 
used. One approach is to set = 30 ~ 1/aw as this is the multiplicity scale at which 
the violation of unitarity arises fiy]. The second approach is to distribute ub according to 



the LOME expression (|^). In either case, the high multiplicity of electroweak bosons means 
that Monte Carlo techniques that can generate energy and momentum configurations of 



many massive particles efficiently have to be used. This is discussed in Sect. 2.1 
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Figure 4: Convolution of step threshold function with proton structure functions for pp 
energies of 40 TeV (dashed line) and 200 TeV (solid line). The parton level cross section, 
(Tt, has been taken as the minimum of either the value from the LOME calculation or the 
S-wave unitarity limit. 



We now give a justification of the assumption that particles can be generated isotrop- 
ically. We saw above that valence quark interactions are dominant. The valence quark 
structure functions at large x have the approximate form 

/ (a;) = (1 - x)" , (18) 



and the convolution of these structure functions with (p!7D produces events with s = X1X2S 
close to the threshold, where higher partial waves are suppressed and so the angular dis- 
tribution is approximately isotropic. 



We can characterize the available energy per particle by a proportionality constant [10 



where 

Mtot = J2^i (20) 

i 

is the total mass of the outgoing particles. The average kinetic energy per particle can then 
be estimated using 

E~(/^-l)m (21) 

for an average particle mass m. For values of fw close to one, the average kinetic energy 
per particle is small, and one can expect that events will be fairly isotropic. This isotropic 
character will be reduced for large values of fw- 
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2.1 Boson final states 



Evaluation of the LOME cross section (|^) has been performed using a Monte Carlo tech- 
nique. This approach is closely related to the generation of events for Monte Carlo sim- 
ulations. Events are generated uniformly in phase space for a fixed particle configuration 
{nB,nH,Cf), where the symbol Cf represents the fermion configuration. For each event, a 
weight that is the product of the matrix element for the event and the phase space integral 
for the configuration is calculated. The average of these weights is then the Monte Carlo 
estimate of the cross section. 

The phase space integral is calculated using the saddle point method . This method 
is based upon Laplace transforming the integral in order to remove the delta function, 
leaving an inverse transform, which can then be evaluated using a saddle point approxima- 
tion. This technique can also be used to evaluate the cross section whilst this allows 
computation of the cross section it is restrictive from an event generation point of view. 

For a given energy and particle configuration, the generation of events is performed 



using the MAMBO Monte Carlo algorithm ||5J] . This algorithm produces configurations in 



phase space obeying energy and momentum constraints. The algorithm is more efficient 



than others, such as RAMBO |55] from which it was derived, because various parameters are 
tuned for each different configuration of particles and total energy. The events generated 
by MAMBO each have an associated weight w. These weights are distributed in the range 
to 1, and so the correct unweighted distribution is obtained by rejecting events with w 
less than a random number generated in this range. 

The evaluation of the matrix element for each event is the third step in the cross 
section calculation. Computation of the factor in (|^) for each gauge boson, averaging over 
the set of generated events, and multiplying by the phase space integral for the particle 
configuration, gives us an estimate of the leading-order cross section {nB,nH,Cf) for 
a given set of final state particles. The total quark level cross section is then obtained by 
summing over all possible particle configurations. The particle numbers are constrained by 
the requirement that the total mass of all the outgoing particles has to be less than the 
energy of the interaction. 

In summing over the number of vector bosons, the relative number of Z" bosons and 
photons has to be considered |Q. The instanton is a field configuration in unbroken SU(2) 
gauge theory; in the standard model this symmetry is broken to give the physical particles 
W^, Z" and 7. Whilst equal numbers of W"*", W~ and bosons of the unbroken symmetry 
are produced by the instanton interaction, the bosons are projected onto the Z° and 
7 physical states, with the relative probabilities of cos^ 9w and sin^ 9w respectively. The 
probability distribution for photons and nz bosons is then described by a binomial 
distribution. 

This procedure also yields the mean boson multiplicity as a function of energy, as 

-UB = (22) 

Ens ^LO {nB, Cf) 

A plot of hb is shown in Fig. ^, along with the estimate from the 'holy grail' function F, 
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Figure 5: Mean number of bosons as a function of energy. Diamonds LOME prediction, 
solid line Eq. (|l|). 



which is shown as a solid line on the same plot. The difference is due to the reduced phase 
space available when one considers massive electroweak bosons. 

At a given multiplicity ns, integration over all possible SU(2) configurations of the 
instanton produces a non-trivial momentum structure in the final state, as noted by Es- 
pinosa 0. This structure contains terms of the form (p.q) and {ef^uapP^iQukal/s), where p, 
q, k and / are particle momenta and polarisation vectors. The large number of particles 
from the instanton interaction means that there will be many (~ [4^^- + 2nB + — 1]!!)0 



such terms |£6[. Whilst calculation of all these terms is not feasible, it is reasonable to 
assume that their overall effect will be to produce an approximately uniform distribution 
in solid angle. However, these terms could modify the energy spectra of the outgoing par- 
ticles. For example, the LOME expression (0) was obtained without performing any SU(2) 
averaging, and in this case the bosons each have a factor of (4i?^ — m^) /m^. This factor 
favours configurations with higher boson energy. Such a factor has not been included in 
the simulations when either the ub = l/cxw or LOME prescriptions were used. 

The averaging over SU(2) configurations also requires that the number of Higgs bosons 
involved in the interaction is even. This can be identified with the observation that the 
integration over all SU(2) configurations of an odd number of rotation matrices is identically 
zero. Recall that the number of Higgs particles uh is expected to be small from (0). We 
model the number of Higgs particles at a given boson multiplicity ns using the relative 
sizes of the LOME estimate (0) function of nu- 



^We use the double factorial notation, whereby n\\ = n{n —2){n — 4). ..3.1 for odd n. 
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2.2 Fermion final states 



The structure of instanton-induced BLNV involves three quarks and a lepton from each 
fermion generation. This structure arises because the instanton determinant includes a 
zero mode for each fermion doublet; there is a quark doublet for each of the three colours 
of QCD. The integral over all SU(2) configurations ensures overall charge conservation, and 
any difference in total charge between the incoming and outgoing fermions is balanced by 
the relative numbers of charged bosons. 

In order to generate the fermion configuration, one member of each doublet has to be 
included in the interaction. The incoming quarks from the protons determine two of the 
three quark doublets for the lightest generation, leaving the other seven antiquark and 
three antilepton doublets to be chosen. 

Our simulations use a simple phase space model to perform the fermion selections. The 
calculation of the phase space for n particles, V^(s), by the saddle point method leads, via 
a Laplace transform, to the expression 

"1 pC-\-'iOQ 

Vn{s) = — / dpp'h (PV^) $ m (23) 

where 

n 

^' (/5) = n <\>^ iP) (24) 

1=1 

is the product of the individual contributions of each particle to the overall phase space, and 
Ji is the entire modified Bessel function. The value of c is chosen so that all singularities 
of the integrand lie to the left of the contour in the complex (3 plane. The contribution of 
a particle of mass rrii is 

where Ki is the singular modified Bessel function, and in the massless case is simply 



(26) 



The integral ( p3D can be performed in the saddle point approximation. Note that the lim- 
iting values of the saddle point, f3, for the cases of all n particles in the extreme relativistic 
and non-relativistic limits are 



ER 



2n - 3/2 



3 - 3 J!i^ .27) 



respectively. 



Neglecting the effect on (3 of choosing between two particles in a doublet with different 
masses, we can estimate the relative probability of picking a particle of mass mi as opposed 
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Figure 6: Relative probability of emission for fermion doublets. Energy = 20 TeV except 
for the lower tb line. 

to one of mass m2 as 

P = 

01 + 02 

\ miKi{mi(jj J 

for two massive particles, and by 

p = (l + m2^Ki (m2^) ) (29) 

for a doublet with mi = 0. 

We find that, with the exception of the tb doublet, the effect of mass on the probability 
p is negligible. For the heavy quark case, the relative probability of t compared to b 
production is reduced as the boson multiplicity increases, as shown in Fig. It is evident 
that the only significant effect occurs at large boson multiplicity, where the phase space is 
sensitive to mass effects. 

The fermion configurations are described in terms of the gauge eigenstates of the elec- 
troweak theory. It is the mass eigenstates of the quark fields that are observed, and the 
gauge eigenstates of the = — ^ quarks have to be projected onto the mass ones. The 
states are related by the CKM matrix U by 

I mass) = f/ 1 gauge) . (30) 
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2.3 Hadronisation of B violating final states 



The objective of a Monte Carlo event generator is to generate complete final states in all 
the detail that could in principle be observed by a perfect detector. In the case of colourless 
produced particles (electroweak and Higgs bosons, and leptons) it is a reasonable approxi- 
mation to suppose that they interact with the detector or decay independently, apart from 
possible spin correlations that have not yet been computed but are probably of secondary 
importance. For coloured particles (quarks and gluons), however, we know that confine- 
ment plays a role that cannot be neglected, giving rise to hadronic jets whose properties 
are only indirectly related to those of the primary partons. One of the main successes of 
modern event generators [^^ -|60|| has been the surprisingly accurate representation of jet 
properties that they give when applied to conventional processes such as e"'"e~ annihilation. 
It is therefore natural to take over the detailed machinery of jet fragmentation in such gen- 
erators and apply it, unchanged as far as possible, to the unconventional process of baryon 
number violation. 

According to the viewpoint adopted in all the widely-used QCD event generators, the 
process of jet production and fragmentation takes place in three distinct phases. First, on 
the very shortest relevant timescale (the electroweak scale in this case), quasi- free partons 
are produced with the distributions prescribed by the hard process matrix elements. Next, 
these primary partons give rise to parton showers via successive gluon bremsstrahlung and, 
less commonly, quark- ant iquark pair production. The showers can be followed perturba- 
tively until the timescale approaches the hadronic scale 1/Aqcd, when the running QCD 
coupling as becomes large. At this stage the shower development is terminated by imposing 
a parton virtuality cutoff Qq > Aqcd- 

The final phase of jet fragmentation, in which the showers at the cutoff scale Qq are 
converted into hadrons, is called hadronization. This process is not well understood at a 
fundamental level but various models have been developed which can describe it in some 



detail. In the HERWIG Monte Carlo program |53, ^ , a cluster model [q^ is used: all 
remaining gluons in the showers are split non-perturbatively into quark-antiquark pairs 
and colourless mesonic clusters form in the resulting cloud of quarks and antiquarks. The 
observed hadrons come from cluster decay according to a simple phase-space model. In 



the JETSET program [q^, ^ a string model |6j is used: instead of clusters, more ex- 
tended colour-singlet objects called strings are formed and these fragment sequentially into 
hadrons. 

In the generation of the parton showers, and also for the formation of clusters or strings 
during hadronization, it is essential to keep track of the colour structure of the process, 
at least to leading order in where N is the number of colours. An effect of special 
importance is the colour coherence of soft gluon emission by the primary partons, which 
sets the limits on parton showering via angular ordering We first review how this 

works for conventional baryon-number conserving processes,^ and then show how it applies 
equally well in the presence of baryon number violation. 

Consider for example the diagram in Fig. representing a contribution to the ma- 



^For a more extensive review, see for example Ref. 
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Fig 7(a) 



Fig 7(b) 



Figure 7: Soft gluon coherence: (a) emission from a qq pair, (b) emission contribution from 
two quarks after a BLNV event. 

trix element squared for emission of a gluon of momentum q from a colour- singlet quark- 
antiquark pair ij. In Feynman gauge (neglecting quark masses) this is the only contribution, 
and in the soft limit it gives a spin- and colour-averaged emission probability of the form 



N Pi-qPj ■ q 2(27r)3g 

= CA^-Bw, (31) 

TT q An 

where t° are the colour matrices in the fundamental representation, Cp = Tr(t"t")/iV = 
(A^^ — 1)/2N and the radiation function Wij is given by 

_ ^_(l^cos%^_ ^ ^ 
(l-cos^,)(l-cos^^,) ' 

dij being the angle between lines i and j and di {9j) the angle of emission of the gluon with 
respect to line i (j). 



The way in which the radiation function ( |32D leads to angular ordering can be seen as 
follows. We write 

W^,=Wi^ + Wij (33) 

where 

W'=^(l + . (34) 



Then W-j and W-j contain the leading collinear singularity as 0^ ^ and 9j respec- 
tively. Furthermore if we average W-j over (pi, the azimuthal angle of emission with respect 
to line i, we find precisely 

That is, after azimuthal averaging, W-j represents emission inside a cone centred on line i 
with half-angle dij. Similarly W-j represents emission in a cone centred on line j with the 
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same half-angle. Outside these cones, the radiation averages to zero, which is what is meant 
by angular ordering. In the 'coherent parton shower' approximation used in Monte Carlo 



simulations ||66|, the radiation function is set equal to zero outside the angular-ordered 
region. This gives the correct results for quantities that are azimuthally averaged, such 
as multiplicity and energy distributions. The approximation is also found in practice to 
be quite accurate for global quantities that are in principle sensitive to azimuths, such as 
event shapes, because the contributions from angular-disordered regions tend to cancel and 
are in any case numerically small, such regions being far from all coUinear singularities. A 
similar analysis can be performed for soft gluon emission from heavy quark systems: here 
the coUinear singularities are screened but angular-ordered regions are still enhanced ||67| . 

The extension of the above treatment to baryon-number violating processes poses no 
special problems since colour, together with other generators of local gauge symmetries, 
remains absolutely conserved in the hard process. We consider the soft gluon emission 
contribution in which the gluon connects lines i and j emerging from a baryon-number 
violating vertex. Fig. |^. For colours, the vertex couples quarks via the invariant 
tensor eijk.... The emission probability becomes 



dP,. = Cb--^W,, (36) 



where the colour factor is 



e^Jk..e^,k... N{N - 1) iV - 1 ' ^ ' 



If we now make the decomposition ( p3D of the radiation function, we see that the full 
contribution with the leading coUinear singularity along the direction of line i is given by 
Eq. ( pID with the replacement 

1 



E • (38) 



N- 



That is, we should simply average the contribution with respect to all other lines emerging 
from (or entering) the same baryon-number violating vertex. In the coherent parton shower 
approximation, this means that we should choose the cone which limits emission from line 
i to have half-angle 6ij where j is chosen at random from all other lines at the same vertex. 
Within this approximation, the parton showers from each line then develop independently. 

As already mentioned, in the cluster hadronization model adopted in the HERWIG 
program all gluons remaining at the end of parton showering are split into quark- 



antiquark pairs. Colour singlet qq combinations are then used to build clusters. This is 
actually done in the leading- 1/A^ approximation, which greatly simplifies matters since each 
external colour line in any diagram is then uniquely connected to an external anticolour line 
with which it forms a singlet. A similar approximation is used in JETSET [^,^] to decide 
which quarks and gluons should be connected by string segments. For baryon-number 
violating processes essentially the same procedure can be adopted: all connected external 
colour-anticolour pairs are first used to form singlet mesonic clusters or strings, leaving sets 
of A^ colour or anticolour lines connected to each other via the same hard vertex, which 
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form singlet baryonic or antibaryonic clusters/strings. Any excess baryons or antibaryons 
naturally appear amongst the decay products of the latter. 

It follows from the above hadronization mechanism that the excess baryons in a B- 
violating final state tend to be amongst the lower -momentum jet fragmentation products, 
rather than the leading particles which carry the bulk of the jet momentum. An excess 
baryon can be a leading particle only if the corresponding excess quark does not emit 
any gluons above the hadronization cutoff Qq. The probability of this is suppressed by a 
Sudakov form factor, which decreases faster than any inverse power of the hard process 
scale. 



2.4 Backgrounds to B and L violating processes 

We consider three possible backgrounds to BLNV processes: non-perturbative multi-W(Z) 
production, multiple electroweak boson production by perturbative standard model pro- 
cesses, and QCD multijet production. 

Non-perturbative multi-W production is expected to occur at high multiplicities, 
l/aw, with an energy scale ~ mw/aw- A detailed discussion, based on the argument that 
at high energies electroweak theory has a large, essentially constant total inelastic cross 
section, may be found in ref. []TU[. The fermion structure for multi-W production without 
BLNV is 

gi + g2 ^ gs + g4 + ub w(z) + uh h, (39) 

giving a signature similar to the BLNV case considered previously (|^), without the primary 
fermions. It is expected that these processes will have a similar threshold behaviour and 
cross section to the BLNV ones. We have therefore modelled this background in the same 
manner as the BLNV signal. In principle, there are also processes initiated by initial-state 
antiquarks, but again, because of the high threshold the qq combination dominates (see 
Fig. For simplicity, we exclude the other incoming combinations, as in the BLNV case. 

Perturbative multiple production of electroweak bosons by standard model processes 
has been estimated in refs. BD|. The dominant mechanism is expected to be multiple 



top quark production: each t decays into a W boson and a 6-quark jet. At 40 TeV energy 
the expected tt cross section is 3 x 10^ fb while that for titi is expected to be about 400 fb. 
In contrast, the predicted direct 4W cross section is only around 6 fb, which is comparable 
to the rate via Higgs bosons decaying into WW. Clearly, in order for the production of 
much larger numbers of W bosons to be observable, a breakdown of perturbation theory is 
required, as assumed above. 

Even though each jet from W decay is in principle well defined, for high W multiplicities 
the number of possible dijet combinations becomes too high for the W^s to be reconstructed. 
At very high multiplicities, the jets will be too close together to be resolved. A typical jet 
cone size, RQ, of 0.7 means that within a rapidity limit of \r]\ < 3 only 127r/(7r0.7^) = 24 
jets can be accommodated without overlapping. Thus, for hadronic W decays, we have to 



^ i? is defined as R — ^/ATf + A(f)^ , -q = —In tan0/2 is the pseudorapidity and is the azimuth, both 
with respect to the beam direction. 
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worry also about pure QCD multi-jet background. 



Perturbative QCD multi-jet production has been studied most in the process of W + n 
jet hadroproduction |]70|,|7l|, since this is an important background for top quark searches. 



The predicted + 2 jet cross section at 40 TeV for jet pt > 50 GeV, ARjj > 0.7 and 
\rij\ < 3 is about 2 x 10^ fb. Tree-graph calculations for 3 and 4 jets suggest a fall in the 
cross section by a factor of 2-3 for each additional jet, as long as phase space limitations 
can be neglected. 

QCD multi-jet production without any direct W boson has been less studied. For the 
same cuts as above, the predicted 2-jet cross section is about 4 x 10^° fb. Four and six jet 
production can be estimated from ref. |]72[ at roughly 4 x 10^ and 4 x 10^ fb, suggesting 



a fall of around a factor of 10 for each additional jet. However, these estimates are based 
on multi-gluon amplitudes, and the W + n iet studies suggest that multi-quark amplitudes 
play a significant role at high multiplicities. This may explain the smaller rate of decrease 
seen in those studies. The Monte Carlo estimates of ref. |73| are indeed somewhat higher, 
but this may be due to a jet definition in terms of invariant mass, rrijj > 6 GeV, which 
corresponds to a very narrow ARjj > 0.1 at = 50 GeV. A fall by a factor between 2 
and 10 per additional jet would imply, for example, a perturbative 30 jet cross section of 
10-9±9 fb. 



In a recent paper [^], approximate expressions for high-order multigluon tree ampli- 



tudes are used to derive the following estimate for the QCD n-jet cross section 

a„~^[^(n,Ar-2 (40) 

where s is the hard scattering energy scale and A = rn'^j^^^^/s. For ln(l/A) ^ n, the 
suppression factor z{n, A) has the form 

z(n,A)~-^ln2(l/A) (41) 



For the jet definitions used above, taking ^/s ~ 10 TeV gives ln(l/A) ^ 10 and z > 1. 
Thus the approach of ref. |]7i| implies a breakdown of QCD perturbation theory. Indeed, 
the instanton approach adopted here for electroweak physics suggests that in QCD such a 



breakdown would be expected for ~ l/a^, with an energy scale of Aqcd/os |]10 



For studies of lepton number violation, the lepton content of the background is relevant. 
Multi-W production and decay would yield many hard charged leptons (about 0.25 per W) 
which would be confused with those from the primary process. In QCD multijet production, 
from Monte Carlo studies of 2-jet events with jet pr > 50 GeV at 40 TeV, the predicted 
number of charged leptons per jet was found to be 0.07, mostly from charm decay. Thus 
30 jets would provide only 2 charged leptons on average, and these would mostly appear 
at low transverse momentum. 

In conclusion, the cross sections for backgrounds to BLNV are as uncertain as that for 
the signal. However, if a signal exists, then B and L conserving multi-W processes are 
likely to occur at a similar rate, and constitute the background most difficult to distinguish 
from the signal. Consequently we shall concentrate on this background in the Monte Carlo 
studies in the next section. 
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Simulations performed 



Energy (TeV) 


hb estimate 


V^o (TeV) 


fw 


E 


17 


1/aw 


5 


1.8 


60 


40 (a) 


1/aw 


18 


7.4 


450 


40 (b) 


LOME 


18 


2.0 


80 


200 


1/aw 


18 


9.3 


580 



Table 1: Monte Carlo simulations performed. All simulations contained a 10^ event sample. 



3 Monte Carlo results 



In this section we describe the results obtained by performing Monte Carlo simulations of 
BLNV processes. The events under consideration are very distinctive, owing to the high 
number of particles involved. In particular, the large number of leptons present and the 
high total transverse energy E^ will be useful experimental criteria for isolating multi-W 
and BLNV events. 

Four cases were studied to investigate the observability of B and/or L violation. The 
configurations of the parameters ub and a/sq used are listed in Table |1], along with approx- 
imate values for the energy fraction fw and the mean energy given by E. This range of fw 
and ub allows useful comparisons to be made. The simulations were performed using an 



event generator, HERB VI |]I2[, which incorporates the parametrizations and assumptions 
described in the previous section. HERBVI is a package designed as an extension to the 
Monte Carlo program HERWIG, which performs the hadronisation of the final state as 



described in Sect. |2.3| . Some further details are given in the appendix. 



The lack of more detailed knowledge of the matrix element means that we have concen- 
trated on the gross properties of the events, such as the excess of antileptons over leptons, 
as opposed to detailed analysis of momentum spectra or rapidity distributions. In each 
of the four cases studied 10'^ simulated events were generated, with m^p = 140 GeV and 
rriHiggs = 300 GeV. 

In order to assess the detectability of BLNV, we have conceived a hypothetical detector 



loosely based on the proposals for detectors at the LHC [[75| ,|76|. The simulations were 
performed with the following acceptance criteria: 

• A minimum of four identified charged leptons of the same family (e or n) for an event 
to be accepted 

• Pseudo-rapidity coverage up to jr/l = 3 

• Minimum particle transverse momentum pt of 10 GeV 

• Electron and muon identification efficiency 95%. 



In addition, we define isolated leptons and photons in the same way as the ATLAS LHC 
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Figure 8: Total transverse energy distribution. 



detector proposal . A particle is considered to be isolated if there is a total of less than 



12 GeV energy deposited in the detector calorimetery around it within a radius R < 0.2. 

No cut on Et has been made in the simulations. In Fig. ^ the distribution of Et for 
the four simulations has been plotted. Here Et is defined as the sum of transverse energy 

1 /2 

(m^ + Pt) for all detected particles. A minimum Et of 5 TeV ( 1 TeV for the 17 TeV 
simulation) will include all the events generated here, yet reject virtually all other types 
of events. Obviously, this cut would have to be 'tuned' in an experimental analysis, but 
combined with the multiple lepton cut the background from any non-multi-W standard 
model process should be essentially zero. 



3.1 Boson Spectra 

With the possible exception of photons, the study of bosons from instanton-induced pro- 
cesses will be difficult. It may be possible to reconstruct the leptonic decays of Z°'s; for 
decays this will not be possible due to the large missing momentum from the neutrinos. 
The identification of jets, and hence the reconstruction of hadronic and Z° and decays, 
will also be difficult because of the large number of overlapping jets. However, the energy 
spectra of the bosons produced in the initial interaction are an important aspect of the 
process since they reflect the characteristic size of the interaction region ~ 1/mw- 

The boson energy spectra for the simulations are plotted in Fig. p. The distributions 
have been weighted so that the total area under each graph is equal to the average number 
of each type of boson produced in the interaction, with the exception of the Higgs boson 
results which have been increased by an factor of ten. The much lower average boson 
energy of the fw ~ 2 simulations compared to the other two can clearly be seen. 

The distributions from these plots are shown again in Fig. |TD| along with the 
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Figure 9: Boson energy spectra. Area under each curve is the average number of bosons 
by type per event. Sohd W^, dotted dashed photons, dot-dashed Higgs (xlO). 
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Figure 10: Comparison of boson energy spectra with sphaleron decay spectrum. 



sphaleron decay spectra |44| for comparison. The sphaleron decay corresponds to fw = 2. 



It can be seen that the distribution differs slightly from that of the fw = 2 simulation. 
The calculation of the sphaleron decay products involves the solving of the non-linear 
field equations. Therefore the spectra of the outgoing particles from this decay include 
the effects of final state corrections (but not of the integration over orientations of the 
classical background field). The similarity of the curves shows that the assumption of a 
flat distribution of outgoing particles in phase space is a reasonable first approximation. 
Note however that this assumption yields more particles in the high momentum tail of the 
distribution. 

Direct photons from the primary interaction, which have not been considered before, 
may provide the best opportunity to investigate the bosonic spectra of the instanton in- 
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Figure 11: Photon energy spectra for the 17 TeV simulation. Sohd hne all photons, dashed 
photons from primary interaction, short-dashed isolated photons. 
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Figure 12: Photon energy spectra for the 40 TeV simulation, using the leading order matrix 
element estimate for ub- Solid line all photons, dashed photons from primary interaction, 
short-dashed isolated photons. 

teraction. In Fig. |Tl] the photon E and pt spectra from the 17 TeV simulation are plotted. 
Typically ~ 60 photons are produced per event, the majority of which are from pion decay 
and at the low end of the energy spectrum. The tails of the distributions largely comprise 
direct photons from the primary interaction. The average number of photons produced in 
this way is sin^ Ow 1^, which corresponds to ~ 2.2 for the ~ 1/aw and ~ 7.7 for the 
LOME estimate respectively. The short-dashed line is the spectrum of photons after the 
isolation cut had been imposed. For this simulation, the domination of the photon spectra 
by primary photons for E > 100 GeV and pt > 100 GeV means that measurement of 
photons in this region would provide information on the BLNV process structure. Similar 
conclusions can be drawn from the other two ri^ ~ 1/ aw cases, but in the 40 TeV simula- 
tion using the LOME estimate, shown in Fig. there are very few isolated photons. 
This is because of the larger number of particles in the detector, which makes finding an 
isolated particle less likely. 

The isolation cut removes the lower energy particles, leaving the tails which are mainly 
composed of photons from the primary interaction. The study of these tails will provide 
useful information about the underlying process. Any structure from the matrix element 
of the instanton process will show up in this part of the spectrum. Therefore measure- 
ments of photon spectra will be important in any experimental study of BLNV or multi-W 
phenomena in general. Good calorimetry, especially if a 'tighter' definition of isolation is 
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Figure 13: Energy and transverse momentum spectra for baryons and antibaryons. Solid 
line, all baryons and antibaryons, dashed antibaryons from BLNV vertex only. 



possible, would aid this task considerably. 



3.2 Bar yon number violation 



The high energy threshold of the BLNV processes we are studying, along with the 'central' 
nature of the events, means that the events are characterised by a large number of particles, 
O ^oo/^~^^, within the detector. Nevertheless, many final state particles escape detection, 
either because they do not enter the solid angle covered by the detector, or because their 
energies are so low that they cannot be reliably detected. Consequently, any verification of 
BLNV will have to concentrate on averaging quantum numbers over a number of events, with 
the aim of demonstrating an asjTumetry. We therefore focus our analysis on determining 
quantities that can be used to demonstrate such asymmetries, and on the number of events 
that will be required. 

The large number of particles present in the detector makes the identification of baryons 
difficult. This is particularly true for any identification technique based upon the recon- 
struction of the masses of unstable baryons. Furthermore, the excess antibaryons from the 



BLNV process are at low energies and transverse momentum, as can be seen in Fig. [13 
The task of separating these antibaryons from the others produced in the interaction is 
effectively impossible. 

We may demonstrate the difficulty of verifying B violation by making an estimate of 
the number of events required for observing B violation with a 95% {3a) confidence level, 
assuming for the moment perfect baryon identification within the detector solid angle. We 
define the baryon number difference Db SiS 

Db = Nb- Nb (42) 

where Nb is the number of detected baryons in a given event, and Nb the corresponding 
number of antibaryons. The averaging is performed separately for BLNV and multi-W 
events. For the 40 TeV /BLNV/ simulation using the LOME estimate for n^, the distribu- 



tions of Nb and Nb are plotted in Fig. 14. The means of these distributions are given in 
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Baryon number results 


Measurement 


BLNV 


Multi-W 


Nb 


15.7±0.04 


15.8±0.04 


Nb 


16.0±0.04 


15.7±0.04 


Db 


-0.26±0.03 


-0.02±0.03 



Table 2: Total number of baryons and antibaryons, within solid angle covered by detector. 
Results from 40 TeV simulation with ub given by the LOME estimate. 
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Figure 14: Number of baryons and antibaryons observed in detector per event. Solid line 
baryons, dotted line antibaryons. 

Table |^. 

In order to estimate the number of events required for baryon number violation to be 
established at a suitable confidence level, we compare two sets of data. The first is a mixture 
of BLNV and multi-W events in equal amounts, corresponding to an equal probability for 
these two processes. The second set of data is the multi-W process alone, and we calculate 
the confidence level of separating the two data sets as a function of the number of events, 
using Student's t test [^], as described in appendix B. 

This analysis predicts that approximately 8.10^ events are needed to verify B violation 
at the 95% confidence level. However, this is assuming a 100% efficiency of baryon identi- 
fication, which is far from realistic. The number of events required for a given confidence 
level is approximately proportional to 1/C^ for an efficiency of (. Using a figure of 1% for 
this efficiency, which is probably a very generous overestimate, the number of events is then 
~ 10^. By concentrating on types of baryons with particularly distinctive decay modes, 
it may be possible to decrease this figure, but it is unlikely that any reduction by orders 

^We thank M.A. Parker for useful discussfons on this pomt. 
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Events passing lepton counting cuts 



Energy 
TeV 


BLNV 


Multi-W 


Electron 


Muon 


Either 


Electron 


Muon 


Either 


17 


55.8 


55.2 


81.1 


38.4 


38.0 


63.2 


40 (a) 


72.5 


70.5 


92.8 


51.1 


49.7 


77.2 


40 (b) 


99.5 


99.6 


100.0 


99.3 


99.3 


100.0 


200 


75.0 


73.0 


94.0 


51.7 


50.9 


77.5 



Table 3: Percentage of events passing lepton count cuts. For the 40 TeV simulations, 
case (a) is ub ~ ^/ctw, case (b) LOME prescription. 

of magnitude is possible. We conclude that the verification of B violation is effectively 
impossible for the types of processes considered here. 



3.3 Lepton number violation signatures 

A more promising signature for BLNV is the violation of lepton number. Unlike the baryons 
from the primary process, which tend to be found amongst the low-momentum debris 



after hadronization of the excess antiquarks as discussed in Sect. |2]^, the excess leptons 
emerge directly from the primary process with high momenta and can be identified. It is 
not possible to demonstrate L violation for a single event, as it is not possible to detect 
the produced neutrinos. However, it may be possible to demonstrate an asymmetry by 
considering the average number of leptons and antileptons produced in a set of events. 
In this section we discuss the possibility of observing lepton number violation using both 
model-independent and model-dependent techniques. We restrict our discussion to and 
/i^ leptons. 

Simulations performed for the four cases listed in Table |l| show that L violation is 
much easier to demonstrate than B violation. Therefore a more detailed analysis of the 
characteristics of the leptons is warranted. Recall that a cut on the minimum number of 
detected leptons has been imposed. In order for an event to be included in the analysis, 
it has to contain at least four electrons or four muons. The cut chosen is similar to those 
that will have to be made at the trigger level at any high luminosity supercollider. This 
cut also helps in some cases to reduce the multi-W background. Table |^ summarises the 
number of events passing this cut. 

The simplest method of verifying the non-conservation of lepton number is to measure 
the quantity 

A = Ni- - Ni+ (43) 

and demonstrate that it is incompatible with lepton number conservation. Such a mea- 
surement can be considered model independent in the sense that it does not depend on the 
details of the matrix element of the hard process, provided the effect of the detector cuts 
at low pt and high rapidity can be neglected. 
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Lepton number results 


Energy 






TeV 


BLNV 


Multi-W 


BLNV 


Multi-W 


17 


-0.48±0.02 


-0.09±0.02 


-0.47±0.02 


-0.07±0.02 


40 (a) 


-0.47±0.02 


-0.09±0.02 


-0.49±0.02 


-0.11±0.02 


40 (b) 


-0.47±0.03 


-0.08±0.03 


-0.45±0.03 


-0.10±0.03 


200 


-0.40±0.02 


-0.04±0.02 


-0.41±0.02 


-0.10±0.02 



Table 4: Average value of lepton number differences De and D^. Results based on 10^ event 
simulations. 
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Figure 15: Total number of observed leptons per event. Solid line BLNV process, dotted L 
conserving multi-W process. 



One antilepton of each family is produced by the instanton interaction (j^). Conse- 
quently, we expect that the average of Di for a large enough number of BLNV events will 
tend to a value of —0.5 for each lepton family. However, there are a number of factors 
which modify this prediction. The branching ratio of a W(Z) boson to a lepton-antilepton 
pair is 11%(3.3%). Hence, for hb electroweak bosons the decay products will include on 
average O.ln^ leptons of each type, with a Poisson distribution. The average incoming 
charge of the hard process is non-zero in both the BLNV and multi-W cases, which results 
in a non-zero average total charge for the outgoing bosons, and consequently a non-zero 
value of Di even in processes which conserve lepton number. In addition, not all the leptons 
will be detected, and there will also be contributions to Di from the decays of heavy quarks 
and antiquarks. The effect of these factors can be seen in the measured values of and 
-D^ listed in Table ^. 

We first consider the general properties of the Ni± and Di distributions before discussing 
the possibility of observing lepton number violation. The results for A*";- + A'";+ for both 



muons and electrons are plotted in Fig. |I5| for each of the four simulations. The difference 
between the BLNV and multi-W cases can be clearly seen and is due to the extra antilepton 
from the primary interaction in the former case. 



The number of leptons produced and detected is primarily determined by n^, with the 
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Figure 16: Difference in lepton number Di for muons and electrons. Solid line BLNV 
process, dotted L conserving multi-W process. 



main secondary effect coming from the detector cuts and inefficiencies. The effect of the 
cuts depends on the 'centrality' of the interaction. If the threshold is large in comparison 
to the total beam energy then the interaction rest frame will tend to be nearly stationary in 
the detector frame. As the beam energy increases relative to the threshold, the interaction 
frame is more likely to be moving in the detector frame and consequently more particles 
are likely to be lost at high rapidity. 



The lepton number differences Di for / = e, /i are plotted in Fig. |T6|. These plots show 
a similar parameter dependence to the Ni plots. The antileptons from the BLNV processes 
cause the differences between the BLNV (solid line) and multi-W (dashed line) plots. This 
difference is also apparent in the average values of the distributions, shown in Table ^. In 
the multi-W case (Df.) and (D^) are close to zero. The small non-zero values of these 
quantities represent the bias from starting with incoming protons, and the other factors as 
discussed above. 

The mean differences (Di) can be used to demonstrate the existence of L violation. In 
the same manner as the baryon case, the number of events required to do this is estimated 
by comparing an equal mixture of BLNV and multi-W events to a set of multi-W events 
alone. The estimates of the number of events required to be inconsistent with the multi-W 
data at the 95% confidence level are given in Table ^ These estimates are written in terms 
of the total number of events (BLNV and multi-W combined); of the order of 10^ events 
are needed to verify the existence of an asymmetry. 

We shall now consider other methods for demonstrating lepton number violation. Our 
aim here is only to investigate techniques that might aid in the verification of lepton number 
violation. In particular, we consider the possibility of enhancing the BLNV event sample 
to improve the determination of Dg and by studying the momentum distribution of the 
leptons, and also other indirect signatures of lepton number violation. We note that any 
analysis based on enhancing the sample using the momenta is dependent on the structure 
of the matrix element, and a quantitative analysis would require more detailed modelling 
of the underlying processes. 



It has been pointed out in [|r^ that tagging events by the sign of the highest energy 
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Number difference events 



Energy/ TeV 




D, 


17 


660 


770 


40 (a) 


890 


670 


40 (b) 


2700 


2570 


200 


900 


1250 



Table 5: Number of events for demonstrating BLNV at a 95% confidence level by measuring 
lepton number differences. Case (a) at 40 TeV is 1/ aw number of bosons, case (b) LOME 
prediction. 

lepton should enhance the asymmetry caused by BLNV. This is because the most energetic 
l~ in any event can only come from the decay of a gauge boson, whereas the highest energy 
could have originated from either boson decay or the primary interaction. We define the 
energy 

Ei+ = Max E (/+) (44) 

for events with 

Max E (/+) > Max E (/") (45) 

together with the corresponding Ei- for events where the most energetic lepton is an 
l~ . Events with no leptons of the type under consideration passing the detector cuts are 
discarded. 

The average energy of a primary lepton from the BLNV interaction, which can only 
contribute to the Ei+ distribution, is ~ whereas a lepton from the decay of a gauge 
boson, which can contribute to both distributions, has an average energy ~ (£■ + TOvy)/2. 

Therefore the enhancement of a sample of BLNV events by selecting according to the highest 
energy lepton will depend on the ratio of E to niw. For multi-W events, all the leptons 
are produced by gauge boson decay and so there is only a small difference between the Ei+ 
and El- spectra due to the non-zero charge of the incoming protons. 

At large n^, it is permissible to make the approximation fh ~ mw This gives 

E^mw fwmw 

—2 2— ^^^^ 

whereas 

Er.mwifw-1). (47) 

In the limiting case of fw close to 1, leptons from the BLNV process will have much less 
energy (close to zero) than those from gauge boson decay (~ mw/'^)- As fw is increased, 
this imbalance is reduced. At fw ~ 2, the average kinetic energy of the BLNV leptons will 
be similar to that of leptons from gauge boson decay. As fw increases above this value, 
the BLNV leptons will tend to have higher average kinetic energies than their boson decay 
product counterparts. Therefore, as fw increases the difference between the Ei+ and Ei- 
spectra will become more pronounced. Note that there will always be some asymmetry 
between the spectra caused by the extra antilepton from the BLNV process. 
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Figure 17: Highest energy lepton spectra for 40 TeV {riw ~ I /aw) simulation. In Figs. 17- 
20, the solid lines represent BLNV events, the dashed lines the L-conserving background. 




Figure 18: Highest energy lepton spectra for 40 TeV (LOME prediction for ub) simulation. 

Using fw to classify the simulations listed in Table |l|, we can divide them into two 
pairs. The distributions of Ei+ and Ei- are plotted for the two 40 TeV simulations in 
Figs. p!7| and p!8| for fw ~ 7.4 (n^ given by the LOME prescription) and fw ^ 2 {ub ^ ^/^w) 
respectively. The enhancement of the Ei+ spectra at large fw can be clearly seen. 



The spectra of Figs. and |I8| have been normalised so that the total area under 
each curve is proportional to the number of events contributing to that distribution. By 
considering the relative entries to the two distributions Ei+ and Ei~ , it may be possible to 
observe an asymmetry between the and /~ data. As shown in Table ^, a small asymmetry 
is expected in the case of L conserving multi-W events, due to the initial state charge, and 
a much larger one in the BLNV case. The estimates for the number of events required to 
demonstrate an asymmetry, using the same mixing of BLNV and multi-W events as before, 
are shown in Table ^. 



28 



Relative numbers of highest energy leptons 



Energy (TeV) Lepton 



17 



40 (a) 



40 (b) 



200 



e' 

e~ 
fx' 
e" 



BLNV events 
58l 
40.4 
59.7 
38.9 
64.8 
34.8 
64.7 
34.8 
55.8 
44.2 
55.9 
44.1 
63.9 
35.6 
63.7 
35.7 



Multi-W events 
50^2 
46.4 
49.7 
47.0 
50.1 
47.7 
50.9 
47.4 
50.6 
49.4 
50.0 
50.0 
49.7 
48.3 
50.1 
48.1 



Table 6: Relative number of entries, as percentages, of the lepton energy distributions 
Ei+ and Ei- for 10^ event simulations. For the 40 TeV simulations, case (a) refers to 
ub ~ l/cKvK) case (b) LOME prescription for ub- 



No. events rel. occurrence 



Energy/ TeV 


/^^ 


e± 


17 


340 


285 


40 (a) 


150 


145 


40 (b) 


1130 


1375 


200 


170 


180 



Table 7: Number of events for 95% confidence level by measuring relative occurrence of 
and as the highest energy lepton of an event, compared with the multi-W background 
predictions. 



29 



Energy asymmetry 


Energy /TeV 


Electron 


Muon 


17 


3.7 10^ 


1.7 10^ 


40 (a) 


6.0 10^ 


9.5 102 


40 (b) 


3.5 10^ 


2.2 10^ 


200 


7.6 10^ 


7.1 102 



Table 8: Number of events for 95% confidence level by measuring asymmetry in energy of 
highest energy lepton by event. 



Number difference events 


Energy/ TeV 


Dt 


D- 






17 


1.1 10^ 


7.5 10^ 


9.4 10^ 


1.1 10^ 


40 (a) 


2.0 10^ 


2.4 10^ 


9.8 10^ 


2.9 10^ 


40 (b) 


1.4 10^ 


1.1 10^ 


1.7 10^ 


2.0 10^ 


200 


1.6 10^ 


2.5 10^ 


4.7 10^ 


4.0 10^ 



Table 9: Number of events for demonstrating BLNV at a 95% confidence level by measuring 
lepton number differences with samples enhanced by the sign of the highest energy lepton. 

The BLNV leptons distort the shape of the Ei+ distributions compared to the Ei- ones. 
Thus, another experimental test is the comparison of these spectra. Within the context 
of our model, we may estimate how many events are needed to see an asymmetry in the 
average values of these spectra, by mixing BLNV and multi-W data in the same way as 
before. The estimates for the number of events are given in Table ^ It should be stressed 
again, however, that the energy spectra are model dependent and so any attempt to verify 
an asymmetry based on analysing the lepton spectra would be less reliable than simple 
lepton counting. 

With this reservation in mind, we now consider refining the lepton counting technique. 
We define two distributions, Z);"^ and D'[ , by separating the measurement of Di according 
to the sign of the highest energy lepton or /^) on an event-by-event basis. These 
distributions are plotted in Fig. |19| for the muons and Fig. ^ for the electrons. The 
corresponding average values of and D'[ for the multi-W case are also non-zero, due 
mainly to the bias introduced by using the lepton charge sign for event selection. The 
averages of these distributions can be used as a signature of L violation, and estimates of 
the number of events required are shown in Table ^ Note that these estimates refer to the 
total number of events. Despite the enhancement of the signal due to the cut on the lepton 
charge, the total number of events required is comparable to working with Di alone, due 
to the fact that we have a reduced sample size. 

Some comments on the experimental considerations of detecting leptons are in order. 
The large number of particles present in the detector as a result of a BLNV process will 



make the identification of the particles difficult. As discussed in Sect. 2.4, there is little 
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Figure 19: Difference in muon number. 




-5 5 -5 5 -5 5 -5 5 
17 TeV 40 TeV (a) 40 TeV (b) 200 TeV 



Figure 20: Difference in eiectron number. 



cfiance of identifying individual jets or reconstructing tfie bosons from tlie underlying hard 
process. However, it may be possible to reconstruct leptonic Z° decays. This would reduce 

the number of Icptons in the Di analyses, and therefore should decrease the size of the 
fluctuations. In addition, the momentum spectrum could then be obtained from the 
reconstruction process. However, a more detailed simulation of the measurement of lepton 
energy and momentum is required to investigate this further. 

The detection criteria used are very simplistic, and any attempt to detect BLNV exper- 
imentally would require a more complete detector simulation. In particular, the isolation 
of electrons is comphcated by the presence of many nearby particles, which may 'fake' an 
electron in the detector. On the other hand, muons are normally identified at the outer 
edge of a detector after having passed through the intervening material. Therefore, muon 
identification will probably be easier than electron identification. 

It may be possible to alleviate some of the detection problems by considering only 
isolated leptons. Repeating the analysis based upon and D^, for the 40 TeV (LOME 
prescription for Ub) simulation leads to an estimate of ^ 10'^ events for a 95% confidence 
level in verifying L violation. The larger number of events required is due to the rejection 
of the non- isolated leptons, which amount to approximately 80%, reducing the statistics 
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considerably. The other three simulations, being based on a ~ approximation, 
have on average less leptons per event; this is offset by the fact that there are fewer particles 
in the detector. In order to investigate isolated leptons further, a more detailed study, 
ideally including a full detector simulation, is required. 



4 Conclusions 

We have discussed the construction, implementation and study of a phenomenological 
model for BLNV processes. The characteristic boson multiplicities and energies which are a 
feature of calculations based upon instanton and sphaleron processes within the Standard 
Model have been incorporated. 

The high boson multiplicities mean that such processes are very distinctive from an ex- 
perimental point of view. Particular characteristics of the processes are their high threshold 
energy, leading to 'central' events in the detector, the spherical nature of the interactions, 
and the high particle multiplicities. The high transverse energy of the final state, combined 
with the large number of high p-r leptons and the inability to resolve jets means that any 
conventional Standard Model background can be easily separated. We have argued that 
the only significant background is that of non-perturbative multi-W(Z) production, with 
typical boson multiplicities ~ 1/aw- 

In our calculations we have taken into account the breaking of the electroweak symmetry. 
We expect that the approximate ratio of outgoing final state bosons in the BLNV process 
is 1 : 1 : cos^Ow '■ srn^Ow for the , W~ , and 7 particles respectively. This has 
allowed the study of the production of photons in BLNV interactions for the first time. The 
high energy and high Pt regions of the photon spectra arc dominated by these photons, 
indicating that the boson distributions from the BLNV matrix elements can be directly 
studied. This study can be enhanced by including reconstruction of leptonic Z° decays if 
possible. 

The studies described here confirm that the demonstration of baryon number violation 
is effectively impossible. The 'extra' antibaryons from the B violating vertex are in general 
produced at low energy and pt, which makes them difficult to isolate from the large number 
of other hadrons produced. In addition, the charge of a baryon or antibaryon is not directly 
related to the baryon number of the particle. There is also the added problem of the low 
efficiency of baryon identification at a supercollider. We quantify our pessimistic conclusion 
by estimating the number of BLNV (and multi-W background) events required to verify B 
violation as at least 10^. 

On the other hand, the studies show that the verification of lepton number violation 
is possible, provided that the cross section of the BLNV process is high enough. Charged 
leptons are relatively easy to identify at a supercollider, and there is the advantage that 
the lepton number of a charged lepton is proportional to its charge. A simple model 
independent test of L violation is to count the number of and l" produced on an event- 
by-evcnt basis, and demonstrate an incompatability with L conservation. We find that 
~ 10'^ events would be required. 
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We have also considered the enhancement of the sample of BLNV events compared 
to the background. In particular, the momentum distribution of the antileptons from 
the BLNV process differs from that of the leptonic decay products of the other outgoing 
particles. Our studies show that using this difference, in general, reduces the number of 
events one needs to demonstrate an asymmetry, to ~ 10^ events. It must be stressed, 
however, that any enhancement of the sample using this requires some knowledge of the 
momentum distribution of the outgoing particles from the BLNV interaction. Therefore, 
any such study is not model independent. 

The predictions made from these studies can be improved in a number of ways. The 
most obvious is to improve the understanding of the momentum distribution and angular 
structure of the outgoing state. Calculation of these quantities for the instanton matrix 
element (or any similar non-perturbative process with a high boson multiplicity) is at 
present not possible due to the large number of terms. Some insight may be gained by the 
type of approach used to study sphaleron decays, which we saw are similar to our model 
predictions in the corresponding region of parameter space. 

We close with some comments on the future prospects for studying not only BLNV 
but non-perturbative processes in general. One possibility is the observation of multi-W 
processes in cosmic ray interactions [|7^, |7U|. In particular, constraints can be placed on 



the cross section for such processes in advance of any supercollider study. Secondly, there 



is theoretical evidence ^^^^ instanton processes in QCD may be observed in deep 



inelastic ep scattering at HERA p2|. 



Appendix A - The HERBVI event generator 



The HERBVI event generator |12| is a package written as an extension to the Monte Carlo 



program HERWIG, and is based on the MAMBO |Q event generation algorithm. Versions 
of HERWIG from 5.7 onwards contain hooks to allow the use of this package. The interface 
to HERWIG consists of four routines: 



HVINIT Initialization routine. This sets up the physical parameters for the simulation. 
Once it has been called, the user can modify these parameters. If it has not been 
called before the first call to HVHBVI, it is executed then. 

HVCBVI Clustering routine. This is executed within HERWIG to enable the antiquarks 
produced in the baryon number violating process to be joined up according to the 
colour constraints described in Sect. ^I3| . 

HVHBVI Main event generator, called by HERWIG for process codes IPROC 7000 to 7999. 
At present, two processes are generated by this routine, 
7000 BLNV event generation 
7100 Multi-W event generation 

HVANAL Initial event analysis. This routine produces a table of all final state particles 
and a flag to indicate if the particle is from the initial BLNV interaction. The rapidity 
and transverse momentum of each final state particle is calculated. 
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In addition, provision has been made for user generation of both the parton-level cross 
section and the boson number distribution. The routines HUSGEN (parton-level cross section) 
and HURBOS (boson number distribution) are called if the appropriate flags (HUFSFG and 
HUFRBN respectively) are set TRUE. 

Final-state particle momenta are generated using the MAMBO algorithm, and to each 
event a weight w is assigned. The HERWIG package produces an unweighted distribution 
by accepting events for which the ratio w/wmax exceeds a random number between and 
1. The maximum weight Wmax is determined by a search at the start of the simulation. 

In order to generate events with a high efficiency, two simplifications are made. Firstly, 
for each event only a relative weight is calculated. This removes the need to calculate 
the phase space integral for a given configuration. This integral is constant for a fixed 
boson multiplicity, and the relative weight for different is taken into account when the 
boson multiplicity distribution is generated. Secondly, events are generated isotropically, 
as discussed above. These two steps reduce the 'tails' of the weight distribution, and 
consequently increase the efficiency of the Monte Carlo. It is straightforward to include 
the effect of a non-trivial matrix element by reweighting the events generated by HERBVI. 
However, this will decrease the efficiency. 

Use of the HERBVI package requires that the dummy routines HVCBVI and HVHBVI be 
deleted from the HERWIG package, and that both programs are compiled with larger com- 
mon blocks. We find that a 10000 event simulation takes approximately 5 to 10 hours to 
run on a Vax4000 workstation. About a third of this time is used by HERBVI when gener- 
ating events. More details of the HERBVI package are available from the authors, who can 
be contacted by electronic mail at gibbsOcbhep . cern . ch and webberQvxcern . cern . ch. 



Appendix B - The Student's t test 

The Student's t test is a standard tool in statistical analysis for estimating the significance 
of a particular set of measurements, or that of the difference between two sets of measure- 
ments. Let xi and X2i with variances al and cr|, be the result of two sets of measurements 
with ni and n2 points in each set respectively. We define the combined variance as 

^2 (d + d\^}l±J^ (48) 
\n2 nij ni+n2-2 

with the number of degrees of freedom 

N = m + n2 - 2. (49) 

The estimate of the significance is then performed by using the Student's t distribution by 
evaluating 

t = (50) 



and computing the confidence level 

V^r(i) 



r(UN 

F{t,N)= I (1 + ^1 dx. (51) 
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F is then an estimate of the statistical significance of the difference between the two sets 
of measurements. A fuller discussion of this test is given in . 
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